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We study bosonic atoms in small optical lattices by exact diagonalization and observe a striking similarity to 
the superfluid to Mott insulator transition in macroscopic systems. The momentum distribution, the formation 
of an energy gap, and the pair correlation function show only a weak size dependence. For noncommensurate 
filling we reveal in deep lattices a mixture of localized and delocalized particles, which is sensitive to lattice 
imperfections. Breaking the lattice symmetry causes a Bose-glass-like behavior. We discuss the nature of 
excited states and orbital effects by using an exact diagonalization technique that includes higher bands. 

PACS numbers: 03.75.Lm, 03.75.Hh 
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I. INTRODUCTION 

Trapped atoms in optical lattices offer a fantastic new sys- 
tem for applications in quantum optics, quantum information 
processing, and as a model system for solid state physics. The 
depth and the shape of the lattice potential, which is propor- 
tional to the square of the laser field, can be controlled with 
great accuracy f\\\. Feshbach resonances allow the tuning of 
the interaction between the neutral atoms in a wide range |2]. 
As a fascinating new development, experiments with a small 
number of lattice sites and in particular finite optical chains 
have become focus of actual research. Prominent examples 
are quantum registers in the context of quantum information 
processingJH 0, Et], the manipulation of single atoms within 
few sites |Q]7and experiments with double well unit cells |01. 
In this context a fundamental question arises: How similar are 
finite systems compared with macroscopic systems? 

In a macroscopic lattice, bosonic atoms undergo the quan- 
tum phase transition from a superfluid phase to a Mott insula- 
tor when the depth of the lattice potential is tuned from shal- 
low to deep. This was first discussed for liquid helium on 
porous media |8], later proposed for neutral repulsively in- 
teracting atoms in optical lattices [9J, and recently observed 
experimentally |[T3|Tll[l2l- Driven by the competition of re- 
pulsive interaction and kinetic energy the atoms localize on 
single lattice sites in the Mott insulator phase. Consequently, 
each site becomes occupied by a fixed integer number of par- 
ticles leading to a crystal-like situation. In the presence of 
disorder the system can also undergo a transition to a Bose 
glass phase iSl [li S [H [H ■ 

In this article, we study the problem of few repulsively in- 
teracting bosonic atoms in finite linear chains and small two- 
dimensional lattices and discuss how finite size effects influ- 
ence the precursors of the Mott insulator and the Bose glass 
transition. The rich physics of the crossover from a double 
well to mesoscopic systems is investigated by exact diagonal- 
ization using a multiband basis, which allows accurate results 
and the discussion of orbital effects which were widely ne- 
glected so far 

We find a surprisingly strong similarity to the localization 
process in macroscopic systems and show that the pair cor- 



relation function is nearly size independent. Moreover, we 
gain an intuitive insight into the excitation spectrum. Finite 
systems offer also an unique possibility to study the effects 
of noncommensurate filling, in which the localization is sup- 
pressed by the equivalence of the lattice sites. We observe 
the formation of an insulating lowest band and in deep lattices 
we find the coexistence of localized and delocalized atoms. 
However, these delocalized states are extremely sensitive to 
lattice imperfections which force the localization of atoms in 
a Bose-glass-like phase. 



II. THEORETICAL MODEL 

The short range interaction potential of ultracold bosonic 
atoms can be approximated by a contact potential g6{r — r') 

with the interaction parameter g — ^^-Cs, where Ug is the 
s-wave scattering length and m the mass of the atoms IitIi . 
Thus, the Hamiltonian including the full two-particle interac- 
tion in a periodic potential Vp is given by 



H = j d\ 7/>t(r) 



2m 



+ Vp(r) + f7At(r)^(r) 



^(r), (1) 



where ?/'(r) is the bosonic field operator We use the poten- 
tial Vp — Vox cos'^{kx) + Voy cos'^ (ky) + Vqz cos^(fcz) that 
is truncated to {rix, 1, 1) sites for chains and {nx,ny, 1) sites 
for two-dimensional lattices. The periodic potential is contin- 
ued at its boundary by a harmonic confinement potential 1 18| 
(see inset of Fig.!!). We model an optical lattice with the pe- 
riodicity a — 515 nm (k — n/a) and vary the depth of the 
lattice from lEp to AOEp given in units of the recoil energy 

^2 1,2 

Ep = -^^^ — 2.16 /ikHz. The quasi-one-dimensional chains 
and two-dimensional lattices have a transversal confinement 
Voy — Vqz ~ 'iOEp and Vbz = 'iOEp, respectively, so that 
the transversal tunneling can be neglected. Exact diagonaliza- 
tion is performed in the Bloch representation of the optical lat- 
tice. By using a few-particle basis the two-particle interaction 
is fully included. The truncation of the basis at a sufficiently 
high energy allows the inclusion of orbital effects [19]. The 
calculations are performed for **^Rb atoms which are present 
in many experimental setups. 
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FIG. 1: (Color online) The lowest coefficients Cn for two (o) and 
eight atoms (+) in a double well potential in dependency on the lat- 
tice depth Vo calculated by exact diagonalization. The lines without 
markers are obtained using the LBA. The inset shows a symmetric 
and an antisymmetric basis wave function. 



III. A PRECURSOR OF THE MOTT INSULATOR IN 
FINITE SYSTEMS 

We start the discussion with atoms in a double well which 
exhibit a crossover reminiscent of the superfluid to Mott insu- 
lator transition. Although the double well with commensurate 
filling represents an intuitive and easy-to-handle model, the 
eigenstates and the spectrum have already a structure similar 
to the studied chains with 3 to 10 sites. A double well with two 
particles can easily be treated analytically with the restriction 
to the lowest band. The Hilbert space separates into states 
with even and odd parity that do not couple. The subspace 
with even parity comprises of the states |2,0) = "^^s^|0) and 

|0, 2) = -^bl^\0), where bl is the creation operator of parti- 
cles in the symmetric state and bis in the antisymmetric state 
(see inset of Fig. [T]!. Since the energy difference between sym- 
metric and antisymmetric single-particle states Xs and Xas is 
twice the tunneling energy t, the energy of the states |2, 0) and 
|0, 2) differs by 2r2 — At+A. The difference in the interaction 
energies A equals g J cfir [xtsi^) ~ Xs using the real space 
representation. The off-diagonal matrix element between both 
states is / = gj d?r Xs (r)Xas(i")- Thus, the ground state is 
given by = cos 6*12,0) + sin 6*10, 2), where 



atan( ). 



(2) 



For very shallow lattices, when t approaches infinity, 6 van- 
ishes and ^'o.vb^o equals |2, 0), i.e., both particles occupy the 
energetically lower symmetric one-particle state. 

In the limit of deep lattices {t — > 0) the difference between 
the symmetric and antisymmetric wave function x^{x) — 
Xas(a^) vanishes (A — > 0) and d approaches — 7r/4. Conse- 
quently, the ground state is given by ^o.vii^oo = "Tf 1^' ^ 
-^|0, 2). Using the creation operators of a particle in the left 



and right well ~ -^(6s ± b\^ the ground state can be 

rewritten as ^o.Vo^oo = ^J^J|0) and it becomes obvious that 
one particle is localized on the left site and one particle on the 
right site. In addition to the localization, a fundamental prop- 
erty of a Mott-insulator-like state is an excitation gap which is 
given here by 2/. The first excited state with an uneven parity 



and the second excited state 5* 



l/2,yo^cx) 



= \{bfTbfm 
are degenerate |20]. Both states are symmetric-antisymmetric 
combinations of doubly occupied sites and thus represent 
particle-hole excitations. In chains and lattices these linear 
combinations of particle-hole excitations build up the excited 
band as discussed further below. 

Some of the coefficients c„ in the expansion of the many- 
body wave function of the ground state are plotted in Fig. [T] 
Already for a double well it is instructive to compare the 
results that are obtained using the lowest band approxima- 
tion (LBA) with numerical multiband calculations. For Vb ^ 
10i?fl the deviations are well noticeable, whereas in deep lat- 
tices the LBA leads to nearly perfect results for filling factor 
V = \. For higher filling factors the total interaction energy 
and consequently the deviations increase, since higher one- 
particle bands are occupied in order to minimize the interac- 
tion energy. Exemplarily, the coefficients of the states |8, 0) 
and |0, 8) for filling factor v — ^ are plotted in Fig. [T] Ad- 
ditionally, the lowest band coefficients |6, 2), |4, 4), and |2, 6) 
contribute. For higher filling factors than v — \ the deviations 
do not vanish in deep lattices, since in that case v interacting 
particles are trapped on each site, leading to a modification 
of the effective single particle orbitals. This has direct impli- 
cations for the Bose-Hubbard model ^ |2ll] which is widely 
used in this context and is commonly restricted to the lowest 
band. Figure [T] shows that for higher filling factors (y > 3) 
this restriction leads to noticeable deviations from a multi- 
band calculation, whereas for low filling factors the deviations 
are small. Experimentally, deviations for higher fillings have 
been observed, e.g., by measuring the on-site interaction en- 
ergy iij 




FIG. 2: (Color online) The momentum distribution of six bosons in 
a quasi-one-dimensional chain with six sites. The crossover from a 
delocalized wave function (Vo < ^Er) to a localized wave function 
(Vc) > lO-Bfl) can be observed. 



3 



80 
70 

60 


" 1] 






1 


— V^Er . 

--- V5Er 

— V^OE^ ■ 

--- V5E^ 


ra 

X 50 

OJ 

7 40 

C\J 

w ou 




I f 

! j 




-| 1' 




u 

i 


20 














10 
n 


fin 







¥ ') 

V '7 

ly 


Hi ill 

V- 4 
V ' / 






10 


20 30 




N -10 




s 








U / t Xo ^ 



^ t d ) — ^ 



10 12 14 16 18 20 



FIG. 3: (Color online) The pair correlation function of six bosons 
in a chain with six sites shows the localization of particles for 
Vb > lO-Efl. Due to the low filling factor v = 1, the deviation from 
results obtained using the LBA (dotted lines) are only noticeable for 
Vo < 5Er. 



FIG. 4: (Color online) The integral over the local correlation function 
for chains with Ns sites and small two-dimensional lattices (6, 9, and 
12 sites) with filling factor v — 1 (in LBA). Inset (a) is a logarithmic 
plot in units of Vo / Er and inset (b) is a double logarithmic plot in 
units of U I J for a chain with ten sites. 



In the following, we extend the double well to a chain with 
Ns — ^ lattice sites and six particles. The momentum distri- 
bution Pp{kx) of the chain is shown for different lattice depths 
Vq in Fig. |2] where ky ~ = 0. For very shallow lattices 
(Vb = lEft) a narrow central peak indicates the delocaliza- 
tion of all particles over the lattice, whereas for deep lattices 
a broad Gaussian momentum distribution is observable. The 
latter can be assigned to particles that are localized in the cen- 
ter of a single lattice site. Although the system is very small 
the similarity to macroscopic experimental results ifTol [l2ll is 
striking. This is a first indication that the localization mecha- 
nism is not very size dependend (see also L24.1 ) which is dis- 
cussed in detail below. 

For Vb — lEfj minor dips in the momentum distribution 
(Fig.|2]l at fc„ = {n/Ns)a* are observable, where a* = 27r/a 
is the reciprocal lattice vector. These dips originate from the 
suppression of standing waves in the confinement with odd 
parity and wavelengths A„ — Ng a/n with n = 1,2, since 
the ground state has an even parity. Increasing the lattice 
depth to Vb = ^En, Bragg peaks located at the reciprocal 
lattice vector appear. At the same time, the central peak drops 
rapidly in height and becomes broader, i.e., the particles begin 
to separate into different wells as the interaction grows rela- 
tive to the tunneling. At Vb = IOE't? the minima are smeared 
out and only a small modulation of momentum density due to 
delocalized particles remains. This progress proceeds with in- 
creasing lattice depth, so that for approximately Vq — 30Eij 
the momentum distribution has a Gaussian shape, correspond- 
ing to completely localized particles [|25ll . 

However, a Gaussian momentum distribution may also 
arise from a superposition of delocalized states and does not 
prove the localization on single lattice sites. Therefore, the 
pair correlation function 



G2(r,r') = (v3t(r)V't(r')V'(r')^(r)> (3) 
p(r') 

at y = z = is studied, which reflects the conditional density 
of — 1 particles if one particle can be found at r'. In Fig. [3] 
the pair correlation is shown for x' ~ — f , i.e., for one particle 
located on the third site. For Vq < 5i?/? other particles can 
be found on the same site, whereas for Vq > lOEji the pair 
correlation vanishes on the third site completely. Hence, that 
site and consequently all sites are occupied with exactly one 
particle in deep lattices. 

The integral over the local correlation function lG2ix,x) — 
^/dx^p(x) JdxG2{x,x) measures the average probabiHty of 
finding two particles at the same position and is consequently 
a good measure for the total spatial correlation of particles. In 
Fig. |4] the local correlation integral is shown for chains with 
Ng sites and filling factor i/ — 1. The calculations are re- 
stricted to the lowest band which is quite accurate due to the 
low filling (see Fig.O. 

Overall we see an exponential decay of the correlation inte- 
gral with increasing lattice depth. In deep lattices the integral 
vanishes which reflects the localization of all particles. For 
Vb < 6Efi the correlation decreases with an exponent —cVq 
[see fit for Ns = 10 in the inset (a) of Fig. HI. Thus, an in- 
crease of the potential barriers causes a relatively strong sep- 
aration of particles in this regime due to the strong overlap of 
wave functions. For Vq > 8Ej^ the particles are located pre- 
dominately on single lattice sites and the extent to neighboring 
sites is small. In this region the exponential decay is weaker 
and the correlation drops with an exponent —c'Vq^'^. It is re- 
markable that our calculations show a crossover between two 
different correlation regimes already for a small sized system. 

In the studied systems we observe a partial loss of spatial 
correlation in the region corresponding to the superfluid phase 



4 



as also found in Ref. fz^. Moreover, we see finite correlations 
above the transition point, which can be understood quanti- 
tatively by using a perturbative ansatz to first order in t/U 
ifm I26I1 . where t is the tunneling matrix element and U the 
on-site interaction energy. The perturbed wave function can 
be written as ^' = ^-^F/ + {f Y.^, b\b,W^jl where ^Ij^ 
is the pure Mott insulator state at infinite lattice depth which 
has a vanishing local correlation. The operator h\hj creates 
a particle-hole state with a hole on site j and a doubly oc- 
cupied site i. For doubly occupied sites the pair correlation 
has a constant value if neglecting the interaction. Hence, the 
expectation value of the correlation integral is roughly pro- 
portional X.Q {t/UY. In a double-logarithmic plot in units of 
U /t [inset (b) of Fig. HJ the correlation integral shows a lin- 
ear behavior between Vq — GEn and 17 En, i.e., the integral 
is proportional to {U/t)^. We observe a value for /3 that is 
slightly above /3 = —2 (about 13%). For shallow lattices the 
simple perturbative ansatz is obviously not suited and for very 
deep lattices (Vq > 20Eji) the expectation values are above 
the power law fit (see also Ref. lUlll '). 

Only in very shallow lattices the correlation integral varies 
noticeably with the number of lattice sites, whereas for Vq > 
3E}i the integral is nearly size independent (even for a double 
well system). Additionally, the differences for Vq < 3Eij are 
quite small for more than four lattice sites. Thus, the localiza- 
tion of particles depends very weakly on the number of lattice 
sites Ns, which indicates that the localization may occur in the 
same manner also in chains of macroscopic size. We conclude 
that the blocking mechanism, which is caused by the tunnel- 
ing prohibiting repulsion, is to a large extend insensitive to 
changes of the system size. Apparently, the coherence length 
in the insulating region drops below the extend of the system, 
which becomes consequently a good representation of a larger 
one. This explains the similarity of the presented momentum 
distributions and experimental results. 

Exemplarily, the energy spectrum of a system with A^^ = 6 
sites, plotted in Fig.|5] shows the formation of narrow many- 
particle bands for deep lattices. The spectrum is in accor- 
dance with the experimentally observed excitation spectrum 
in Ref. llT2ll . The bands are gapped by the interaction energy 
U of two particles on the same site, so that, for example, the 
first excited band consists of states where one particle inter- 
acts on average with one other particle on the same site. The 
number of states in each band is given by the possibilities to 
remove a certain number of particles and put them onto other 
sites resulting in an interaction energy of nil (e.g., the first 
band has 30 states due to 6 possible sites for a hole and 5 pos- 
sible sites with double occupation). In particular, the energy 
of the nondegenerate ground state decreases and the ground 
state becomes separated from the excited states which is char- 
acteristic for an incompressible Mott insulating state. In the 
limit of deep lattices the interaction energy of the ground state 
vanishes, since the wave function overlap decreases due to the 
localization on single sites. The energy of the excited bands 
slowly increases due to the stronger confinement. Since the 
sites are equivalent, the eigenstates of the bands are delocal- 
ized. These delocalized states that form the excited bands in 
this commensurate system reappear within the ground-state 
band of noncommensurate systems, which are discussed in 



the next section. 

The results obtained so far for momentum distributions, 
correlation functions, and energy spectra of quasi-one- 
dimensional chains can be generalized to lattices. As an ex- 
ample, we study a quasi-two-dimensional lattice with 3x3 
lattice sites with integer filling factor v = 1. The momen- 
tum distribution is presented in Fig. |6]for different potential 
depths. It shows the crossover from delocalization to local- 
ization and compares well with experimental results ll26ll . For 
^Er and 7.5Er Bragg peaks at k = ±a*k^ and k = ±a*ky 
appear in the distribution. Due to the finite number of sites per 
dimension, additional dips at A; = ^a* can be observed. In- 
creasing the lattice depth, the momentum distribution smears 
out due to the localization of particles. This starts at approx- 
imately IQEj^ and is far advanced at 12.5i?fl, which is in 
accordance with the critical point for infinite systems iEtIi . 
The two-dimensionality of the system becomes apparent for 
Vq between lOEfj and 15Efj. The momentum distribution 
reflects the square symmetry of the boundary of a single lat- 
tice site. At roughly Vq — iOEn the distribution becomes a 
Gaussian which indicates that the particles are located deep 
in the wells. Then the confining potential experienced by the 
particles is circular symmetric and can be approximated by a 
two-dimensional harmonic oscillator, which is reflected in the 
momentum distribution. 

A more detailed picture of the localization process is given 
by the integral over the local correlation function la^ix.x) 
which is plotted in Fig. H) Compared with the quasi-one- 
dimensional chain, the localization process is shifted towards 
deeper potentials because of two reasons: The interaction en- 
ergy on each lattice site is diminished due to smaller confine- 
ment in the y direction and the tunneling is enhanced since 
tunneling to four nearest neighbors is possible. For the two- 
dimensional lattice, we observe the formation of gapped ex- 
cited bands in deep lattices as well as the separation of the 
ground state in the energy spectrum in the same way as dis- 
cussed for chains. 




FIG. 5: The energy spectrum of six bosons in a chain with six sites 
relative to the energy Eq of the noninteracting system. The spectrum 
shows the ground state and the two lowest bands. 
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FIG. 6: (Color online) The momentum distribution (normalized to 1) 
of a two-dimensional 3x3 lattice with nine atoms for various lattice 
depths Vo using the LBA. 




We are aware that the finite systems discussed here cannot 
exhibit macroscopic phases and a phase transition. Nonethe- 
less, the precursors to the superfluid phase and the Mott insu- 
lator map on many aspects known from macroscopic systems. 
The momentum distribution (Figs. |2] and |6l), the coiTelation 
function (Figs.[3]and|4]i and the energy spectrum (Fig.|5]) show 
that the localization process is in good agreement with infinite 
systems. Finite size effects can be observed but do not dom- 
inate the behavior of the system. Therefore, simulations with 
few lattice sites are also quite applicable to larger systems and 
offer an intuitive and detailed insight due to the accuracy and 
the inclusion of orbital effects. In the next section we use our 
method to examine systems with noncommensurate filling. 



IV. NONCOMMENSURATE FILLING 

When a noncommensurate filling of the lattice is present 
the physical situation becomes more complicated. The local- 
ization of all particles as in the Mott-insulator-Iike regime is 
suppressed by the symmetry of the potential, since the equiv- 
alence of sites requires the same filling on all sites. Conse- 
quently, particles which in principle would prefer localiza- 
tion must delocalize over the whole lattice. The differences 
to commensurate filling can be illuminated by considering a 
double well with three atoms restricted to the lowest band. 
The basis of even parity states comprises of the two states 
|3,0) = ^^1^|0) and |1,2) = ^6lbl?|0). The difference 
in energy between both basis elements is 2J7 = 4t + A with 
A = gjd^r [xiir) + ixUr)xi{r) - 3x^(r)] and the off- 
diagonal matrix element is / = \/3g f (Pr Xs (r)Xas(i")- Thus, 
the solution for the ground state is given by the same expres- 
sion as for two particles, i.e., vJ/q = cosf?|3, 0) + sin^jl, 2). 
In deep lattices the tunneling energy t vanishes, but A, which 
had vanished for two particles, becomes 4// VS. 

The ground state in the limit of deep lattices is given by 
%yo-oo - ^|3,0> - i|l,2> = ^blblibl +bl)\0), since 
6 becomes 2tt/3. This represents a wave function with two 
atoms that are localized and one atom that is delocalized be- 
tween both wells. Consequently, the ground state is a mixture 
of localized and delocalized particles. Of great importance in 
this context is that in this limit the first excited state, which 



FIG. 7: (Color online) The momentum distribution of TV = 5 to 10 
atoms in a chain with five sites at Vo ~ iOEn. For a perfectly sym- 
metric potential the localization of all particles can only be observed 
at integer filling (dotted lines). 



has odd parity and is givenby = ^bjbKbj — bl)\Q), 

becomes degenerate with the ground state. Small asymme- 
tries lead to linear combinations of the two quasidegenerate 
states and result in the nonsymmetric states -^5|^6j|0) and 

-^bjbl'^\0). Consequently, in a potential with broken sym- 
metry the third particle can localize in one of the wells, if the 
lattice is deep enough. This localization process depends on 
the potential difference between the two wells compared to 
the energy difference between ground and first excited state 
Ei{Vo) — Eo{Vo). The third and fourth excited state which 
form the "first excited band" are separated from the ground 
state by 4//%/3. 

These intuitive results for a double well transfer nicely 
to chains and even two-dimensional lattices. Exemplarily, a 
chain with Ng = 5 sites filled by 5 to 10 particles is studied. 
The momentum distribution for deep lattices (Vo = 40i?fl) 
obtained by exact diagonalization is shown in Fig. [T] For 
integer filling factors (N = 5 and N — 10) a Mott insu- 
lator momentum distribution can be observed as discussed 
in the previous section. When, for example, adding a sixth 
particle to a chain with filling factor = 1, the additional 
particle cannot localize, since all lattice sites are equivalent. 
Despite being delocalized at different sites, the particle has a 
high probability density at the lattice site centers. Therefore, 
the momentum distribution reflects the lattice structure very 
clearly. Correspondingly, it shows peaks at and ±a* as well 
as smaller peaks at ±na* with n — 2,3, which originate 
from the delocalized particles. Additional to this peak struc- 
ture the momentum distribution has an underlying Gaussian 
background, which arises from localized particles. Recapit- 
ulating the ground state for three particles in a double well 

*o,yo^oo = ^^1^1(^1 + ^J)|0)' the interpretation is straight 
forward: For noninteger filling factors ly, the number of par- 
ticles may be written as = kNs + N^^dd with the corre- 
sponding integer filling factor k and the number of additional 
particles A'add- In deep lattices kNs particles localize in the 
wells of the lattice and the remaining A^add are delocalized. 
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FIG. 8: The energy spectrum for six particles on five sites shows the 
formation of a lowest many-particle band. 



The plotted momentum distribution for Vq = AOEji in 
Fig. |2] shows that for = 6 particles the Gaussian back- 
ground is noticeably smaller than for five particles. This in- 
dicates that the localized particles are influenced by the de- 
localized hopping particle which experiences the same repul- 
sive interaction on all lattice sites and consequently interacts 
with all localized particles. For seven particles the height of 
the peaks increases due to two hopping particles, whereas for 
eight particles the background increases, since more sites are 
doubly occupied. Finally, for nine particles all sites except 
one are doubly occupied which is equivalent to the tunneling 
of a hole in a lattice with filling i/ = 2. We also find the delo- 
calization for noncommensurate filling by analyzing the pair 
correlation function 6*2(2;, x') [28]. 

The energy spectrum for six particles in a chain with Ng = 
5 sites is shown in Fig. [8] For deep lattices the formation of 
bands can be observed, which are roughly separated by the 
two-particle interaction energy U. An interesting feature is 
the splitting of the second excited band which would not be 
observable in the LBA. The states with higher energy have 
three doubly occupied sites whereas the other states have one 
triply occupied site. Their energy is reduced by a stronger 
deformation of the wave function accounting for the effec- 
tive repulsive potential created by the two other atoms at the 
same site. As seen before for a double well, the ground state 
becomes quasidegenerate and lies within a band of Ng states. 
Therefore, the ground state is extremely sensitive to small per- 
turbations of the lattice potential, which is discussed in the 
next section. 

It is hardly suprising that the rich physics of noncommen- 
surate filling can also be found in two-dimensional lattices, 
including the formation of a lowest band with its discussed 
implications. For a two-dimensional lattice it is interesting to 
explore the quantum mechanical nature of the states that are 
contributing to the lowest band. For a 3 x 3 lattice with ten 
particles the lowest band consists of nine states due to the one 
additional particle in comparison with commensurate filling. 
Similar to Fig. [8] these states become quasi-degenerate with 
increasing lattice depth representing all delocalized combina- 
tions with the two-particle interaction energy U. In Fig. |9] the 



kx'- , „ , , „ , , „ , , „ , 

FIG. 9: (Color online) The momentum distribution (normalized to 1) 
of the lowest band states of a two-dimensional 3x3 lattice with ten 
atoms at Vb = 30Er (in LBA). The seventh excited state, which 
corresponds to the sixth excited state (rotated by n/2), is not shown 
here [29]. 



momentum distribution of the lowest band states is plotted at 
Vq = 30Eji (the seventh state is not shown) ||29|]. The strong 
interference pattern of the momentum distribution is in eye- 
catching contrast with the Gaussian shape of a Mott insulator 
state and reflects the delocalization. The general structure of 
the shown band states with only few states can be generalized 
to larger noncommensurately filled lattices and to particle- 
hole excitations for commensurate filling. 



V. NONCOMMENSURATE FILLING EM A HARMONIC 
CONFINEMENT 

As a last point we discuss the noncommensurate filling un- 
der the influence of an additional symmetry-breaking poten- 
tial. Perturbations of the potential can cause linear combi- 
nations of the lowest band states that allow the localization 
of the additional particles on specific sites. Consequently, a 
small external confinement can destroy the partly delocalized 
phase in the same way as small random site offsets caused 
by lattice fluctuations in a Bose glass |8]. The parameter 
which triggers the localization of the additional N^^i particles 
is the bandwidth of the lowest band which must be similar or 
smaller than the site offsets to attain localized particles. In ex- 
perimental setups perfectly flat potentials are hard to achieve, 
due to the finite waist of the laser beams, which establish the 
optical lattice, and additional external fields. We investigate 
this effect by using a chain with five sites and six particles 
that experience an additional potential — 2VTe~ '"'o with 
Vt — 40i?i? . The potential is motivated by the Gaussian beam 
waist Wo given in units of the lattice constant a, but it can 
also be approximated by the harmonic potential ^-^x^. The 
momentum distribution at Vq = 20Eji is shown in Fig. [TO] 
for different beam widths wq ranging from wo = 10a to 
Wo = 20Ga. The coiTesponding offset energies relative to 



the central site are given by ei = 



1/ ' 



jE]^ and 62 — 4ei 



(see inset of Fig. [TOb . The width of the ground state band at 
Vb = 20£:fl is roughly ^Er. 

In order to localize the additional particle, the offset ei is 
important with respect to the bandwidth. For wq = 200a and 
Wo — 120a peaks due to delocalization can be well identified 
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FIG. 10: (Color online) The momentum distribution of six atoms in 
a chain with five sites at Vo = 20Eii for a laser beam with a finite 
Gauss width wo. The inset shows the corresponding site offsets ei 
and £2. 



in the momentum distribution, but at wq = 80a (ei — jqEji) 
this structure is smeared out. At wq = 40a (ei — jqEj^) 
the momentum distribution matches the distribution for com- 
mensurate fining shown in Fig.|2](at Vq « 15Eji), i.e., the 
sixth particle is locaUzed in the center of the lattice. Be- 
low Wq = 40a the momentum distribution does not change 
noticeably. Instead, the density changes drastically, since 
the gain in potential energy €2 exceeds the repulsive inter- 
action of two particles at the same site U. At wq = 20a 
the density at the two outer sites vanishes and the inner sites 
are doubly occupied, due to a total energy gain of roughly 
2(62 - [/) ~ 2 X {I.6E11 - 0.6Er). For the stronger con- 
finement Wo — 10a the central site is even occupied by four 
particles. 

Despite the harmonic potential the lattice symmetry can be 
broken by applying different moderate perturbation potentials 
(such as single site offsets or a linear potential). The localiza- 
tion process remains qualitatively the same for the studied fi- 
nite systems, since basicly local site offsets are responsible for 
the localization. We conclude that in perturbed lattices a Bose- 
glass-like localization of all particles can always be achieved 
in sufficiently deep lattices triggered by the ratio of offset 
energies and bandwidth. However, in intermediate lattices a 
mixed phase with A^add delocalized particles can be observed, 
if the lattice fluctuations are smaller than the bandwidth. In- 
creasing the lattice depth (decreasing the bandwidth) first kNs 
particles and in deeper potentials the remaining A^add particles 
localize. Experimentally the observation of the mixed phase 



may be hindered by the finite temperature of the BEC. We 
note that in stronger harmonic confinements in which the off- 
set energies match the interaction energy U (wq < 20a for the 
system above) a precursor of a shell structure can be observed, 
which shows regions with different occupations per site Ii30il . 

VI. CONCLUSIONS 

We have studied bosonic atoms in finite optical lattices us- 
ing an exact treatment which includes effects of higher bands. 
Due to the equivalence of sites, finite lattices with integer fill- 
ing factors exhibit a fundamentally different behavior than 
those with noninteger filling factors. The well studied su- 
perfluid to Mott insulator transition can be recovered in fi- 
nite commensurately filled chains and two-dimensional lat- 
tices with few lattice sites. Our results show the localiza- 
tion of atoms in deep lattices and reveal a striking similarity 
to the momentum distribution observed in macroscopic sys- 
tems. Furthermore, we have shown that the local correlation 
is widely independent of the system size which indicates that 
the localization process in small systems compares with that 
in infinite systems. The localization is also reflected by the 
formation of an energy gap in the energy spectrum. 

For noninteger filling factors only the particles that coiTe- 
spond to integer filling localize in deep lattices whereas the 
additional particles are delocalized. The coexistence of lo- 
calized and delocalized particles in the ground state can be 
observed in the momentum distribution and the pair correla- 
tion function. The energy spectrum shows the formation of a 
naiTow lowest band in deep lattices causing the ground state 
to be extremely sensitive to perturbations of the potential such 
as lattice imperfections or additional confinements. Triggered 
by the ratio of bandwidth and site offsets one observes the lo- 
calization of all particles which is similar to the localization 
process in a Bose glass. In weakly confined systems this leads 
to a localization which occurs in deeper potentials than in lat- 
tices with commensurate filling. 

Briefly, we have shown how the macroscopic physics of the 
Mott insulator and of the Bose glass transfer to finite systems 
and that the detailed simulation of small systems offer impor- 
tant information about larger ones. Especially, we gained de- 
tailed insight into the localization behavior of experimentally 
relevant finite systems. 



VII. ACKNOWLEDGMENTS 



We thank Frank Deuretzbacher for valuable discussions. 



[1] P. lessen and I. Deutsch, Adv. Atom., MoL, Opt. Phys. 37, 95 
(1996). 

[2] S. Inouye, M. R. Andrews, J. Stenger, H.-J. Miesner, D. M. 
Stamper- Kurn, and W. Ketterle, Nature (London) 392, 151 
(1998). 

[3] D. Jaksch, H.-J. Briegel, J. I. Cirac, C. W. Gardiner, and 
P. Zoller, Phys. Rev. Lett. 82, 1975 (1999). 



[4] G. K. Brennen, C. M. Caves, P. S. Jessen, and I. H. Deutsch, 
Phys. Rev. Lett. 82, 1060 (1999). 

[5] D. Schrader, I. Dotsenko, M. Khudaverdyan, Y. Miroshny- 
chenko, A. Rauschenbeutel, and D. Meschede, Phys. Rev. Lett. 
93, 150501 (2004). 

[6] Y. Miroshnychenko, W. Alt, I. Dotsenko, L. Forster, M. Khu- 
daverdyan, D. Meschede, D. Schrader, and A. Rauschenbeutel, 



8 



Nature (London) 442, 151 (2006). 
[7] J. Sebby-Strabley, M. Anderlini, P. S. lessen, and J. V. Porto, 

Phys. Rev. A 73, 033605 (2006). 
[8] M. P. A. Fisher, P B. Weichman, G. Grinstein, and D. S. Fisher, 

Phys. Rev. B 40, 546 (1989). 
[9] D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and R Zoller, 

Phys. Rev. Lett. 81, 3108 (1998). 
[10] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, and 

I. Bloch, Nature (London) 415, 39 (2002). 
[11] F. Gerbier, A. Widera, S. Foiling, O. Mandel, T. Gericke, and 

1. Bloch, Phys. Rev. Lett. 95, 050404 (2005). 
[12] T. Stoferle, H. Moritz, C. Schori, M. Kohl, and T. Esslinger, 

Phys. Rev. Lett. 92, 130403 (2004). 
[13] B. Damski, J. Zakrzewski, L. Santos, P. Zoller, and M. Lewen- 

stein, Phys. Rev. Lett. 91, 080403 (2003). 
[14] L. Fallani, J. E. Lye, V. Guarrera, C. Fort, and M. Inguscio, 

Phys. Rev. Lett. 98, 130404 (2007). 
[15] R. Pugatch, N. Bar-gill, N. Katz, E. Rowen, and N. Davidson, 

e-print arXiv:cond-mat/0603571 (2006). 
[16] R. T. Scalettar, G. G. Batrouni, and G. T. Zimanyi, Phys. Rev. 

Lett. 66,3144 (1991). 
[17] A. J. Leggett, Rev. Mod. Phys. 73, 307 (2001). 
[18] The confining potential V(x) — 7V01 (2^) continues the pe- 
riodic potential at |x| = xo, where n is the number of sites. 

The potential is continuously differentiable which determines 7 

(~ 1) and xo (~ na/2). 
[19] We used up to 100 000 basis states and parity conservation if 

applicable including the formed Bloch bands and the lowest 

"continuum" states (bound states in the confinement). 
[20] The first excited state 'I'l = |1,1) = btbt\0) remains 

unchanged and the second excited state is given by \£'2 = 



cosei'|2,0) +sin6''|0, 2) with 6' = atan[{n + V^V~fP) / 1]. 

[21] D. Jaksch and P Zoller, Ann. Phys. (N.Y.) 315, 52 (2005). 

[22] We can confirm the experimental observation and the calcula- 
tion in Ref. |23] also by the diagonalization of several bosons 
on a single site. The on-site energy U of filling factor v — 5 is, 
e.g., roughly 17% smaller than for = 1 at Vo = 30Er. 

[23] G. K. Campbell, J. Mun, M. Boyd, P Medley, A. E. Leanhardt, 
L. G. Marcassa, D. E. Pritchard, and W. Ketterle, Science 313, 
649 (2006). 

[24] R. Roth and K. Burnett, Phys. Rev. A 67, 031602(R) (2003). 

[25] The density on one site corresponds to the density of a single 
particle in a single-sited cos^ potential, but deviates slightly 
from the harmonic approach. 

[26] I. B. Spielman, W. D. Phillips, and J. V. Porto, Phys. Rev. Lett. 
98, 080404 (2007). 

[27] N. Elstner and H. Monien, Phys. Rev. B 59, 12184 (1999). 

[28] Additional information can be obtained from the pair correla- 
tion function G2{x, x') in deep lattices with one particle fixed 
in the middle of a specific site. As seen before, the correlation 
vanishes at that site for filling factor v = 1. For noninteger fill- 
ing (6 < A'^ < 9) the height of the correlation function varies 
for each site. Thus, regarding the integer number of particles, at 
least some particles must be delocalized between the sites. 

[29] The first and second excited state, the third and fourth excited 
state as well as the sixth and seventh excited state are exactly 
degenerate. 

[30] G. G. Batrouni, V. Rousseau, R. T. Scalettar, M. Rigol, A. Mu- 
ramatsu, P. J. H. Denteneer, and M. Troyer, Phys. Rev. Lett. 89, 
117203 (2002). 



